Orbital Structure and Magnetic Ordering in Layered Manganites: 
Universal Correlation and Its Mechanism 



O 

o 
o 

(N 



I 



I 

o 



O 
O 

o 



I 

o 
o 



X 



S. Okamoto, S. Ishihara, and S. Maekawa 
Institute for Materials Research, Tohoku University, Sendai 980-8577, Japan 

(February 1, 2008) 

Correlation between orbital structure and magnetic ordering in bilayered manganites is examined. 
A level separation between the 3d^,,2_j.2 and 3dj.2_j,2 orbitals in a Mn ion is calculated in the ionic 
model for a large number of the compounds. It is found that the relative stability of the orbitals 
dominates the magnetic transition temperatures as well as the magnetic structures. A mechanism of 
the correlation between orbital and magnetism is investigated based on the theoretical model with 
the two eg orbitals under strong electron correlation. 

PACS numbers: 75.30.Vn, 75.30.Kz, 71.10.-w, 75.80. +q 



Since the discovery of the colossal magnetoresistance 
(CMR), studies of manganites with cubic perovskite 
structure have been renewed theoretically and experi- 
mentally. Competition and cooperation between spin, 
charge and orbital degrees of freedom as well as lat- 
tice cause the dramatic changes of transport and mag- 
netic properties. Manganites with bilayered structure 
^2-2a;5i+2KMn207, whcrc A and B are trivalent and 
divalent cations, respectively, are another class of CMR 
materials . Since an extremely large MR is observed 
near the transition from paramagnetic (PM) insulator to 
ferromagnetic (FM) metal, it has been considered that 
several concepts proposed in the cubic compounds are 
applicable to the bilayered ones. 

In cubic manganites, one of the key factors dominat- 
ing the magnetic orderings is the tolerance factor j^]; a 
bending of a Mn-O-Mn bond decreases the hopping in- 
tegral of carriers. As a result, the ferromagnetic tran- 
sition temperature Tc decreases in the double exchange 
(DE) scenario. However, in bilayered manganites, the 
Mn-O-Mn bond angle is almost unchanged with chang- 
ing cations and carrier concentration, as shown later, in 
spite of a wide variety of the magnetic structures. Various 
key factors dominating the magnetic ordering, which are 
not included in the DE model, were experimentally sug- 
gested, e.g. the antiferromagnetic (AFM) superexchange 
(SE) interaction [Q , the local lattice distortion [||-D , the 
charge and orbital degrees of freedom and their orderings 
P,p^ and so on. However, systematics in their correla- 
tions for a variety of compounds and their mechanisms 
still remain to be clarified. 

In this letter, we study the correlation between mag- 
netic ordering and orbital structure in bilayered mangan- 
ites. The two eg orbitals, i.e. the Sd^z^-r^ and Sd^^-y^ 
orbitals in a Mn'^+ ion split in the crystalline field of the 
bilayered structure and one of them is occupied by an 
electron. It is known that the occupied orbital controls 
the anisotropy of the magnetic interaction as well as its 
strength. The level separation between the orbitals is 
calculated in the ionic model for a large number of the 
compounds. We find a universal correlation between the 



relative stability of the orbitals and the magnetic tran- 
sition temperatures as well as the magnetic structures. 
A mechanism of the correlation is investigated based on 
the theoretical model with the Cg orbitals under strong 
electron correlation. 

We first show that neither the tolerance factor nor the 
bond length governs Tc and the Neel temperature T/v 
for the A-type AFM ordering [^. The tolerance factor 
in the bilayered crystals is defined by t = ((io(i)-A(i) + 

do(2)-A(2))/(2V2dMn-0(3)) where dA-B is a bond length 
between A and B ions. The position of each ion is shown 
in the inset of Fig. |l|(a). Being based on the structural 
data obtained by the neutron and x-ray diffraction ex- 
periments fl^-j2^ , we evaluate t and the bond length be- 
tween nearest neighboring (NN) Mn ions in the ab plane 
'^Mn-Mn ^'^^ ^ Variety of compounds. Tc and Tn are 
plotted as functions of t and ^Mn-Mn Figs. 1 (a) and 
(b), respectively. Almost all t's are located in a narrow 
region where Tc's and Tn's are distributed randomly. 
In addition, Tc is not correlated with d^^_j^^, either. 
Although T/v increases with increasing ^Mn-Mn' ^^^^ 
relation is opposite to that predicted by the DE scenario 
where, with increasing the bond length, the hopping in- 
tegral decreases and the FM interaction in the ab plane 
decreases. We conclude that the DE model, which in- 
cludes the change of the hopping integral caused by the 
change of the bond angle/length, can not explain Tc[n)- 
We also examined correlations between Tc( n) and a num- 
ber of other quantities: the torelance factor evaluated by 
the ionic radius, a Mn-0(3)-Mn bond angle, a Mn-O(l)- 
Mn bond length, Mn-0 bond lengths, a lattice spacing 
between NN bilayers, lattice constants and the valence- 
bond sum for a Mn ion. However, there are not clear 
correlations between these parameters and Tc(n)- 



Let us focus on the correlation between T, 



C(N) 



and 



a relative stability of the Cg orbitals. We employ the 
ionic model to examine the electronic energy-level struc- 
tures. This model may be justified by the following 
considerations |2j]: (1) the manganites at x = are 
classified as charge-transfer type insulators in which the 
ionic model provides a good starting point (2) The 
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FIG. 1. Tc and Tjv as functions of (a) t and 
d^n-Mn- Filled and open circles indicate Tc and Tjv, 
spectively. t and d^n-Mn obtained from the structu 
data in the following compounds: A: Lai,4Sri,6Mn207 
[ii, B: Lai.aSri.TMnaOT (p) Jll, C: Lai.aSri.gMnaOr 
[12, D: Lai.iSri.9Mn207 is) E: Lai.o4Sri.96Mn207 

[[I, F; LaSraMnaOr (s) |i|L_G: Lai.4Sri.6Mn207 (p) ||, 
Pri.4Cai.3Bao.3Mn207 (pIM, I: Ndi.4Cai.6Mn207 {p)_f] 
J: Lai.4Sri.6Mn207 (pIM, K: Lai.2Sri.8Mn207 (p) 
Lai.2Sri.8Mn207 (s) jl^rM: Lai.2(Sro.8 Cao.2)i.8Mn207 
0, N: Lai.2(Sro.7 Cao.3)i.8Mn207 (s) 0, O: Lai.2(Si 
Cao.4)i.8Mn207 (p) [0, P: (Lao.8Ndo.2)iL2Sri.8Mn207 
0, Q: (Lao.6Ndo.4)i.2Sri.8Mn207 (s) @, R: Smi.2Si 
Mn207 (p) S: Lai.2Sri.4Cao.4Mn207 (p) jTl, 

NdSr2Mn207 (p) 0, U: Ndi.iSri.9Mn207 (p) pol, 
LaSr2Mn207 (p) W: LaSri.6Cao.4Mn207 (p) 

NdSr2 Mn207 (s) |22[| where (s) and (p) indicate the s 
gle and polycrystalline samples, respectively. The inset of 
shows a schematic picture of the bilayered structure. 



ionic property is predominant between bilayers. (3) 1 
energy-level structure given by band-structure calcu 
tions shows the same tendency with those by the ionic 
model The energy levels of the Cg orbitals split 

due to the electrostatic potential and one of the orbitals 
is occupied by an electron in a Mn"^+ ion. 

By using a large number of the structural data [p"2|-p3[ , 
we calculate the Madelung potential for a hole in the 
3^322 _j.2 and 'idx'^-y^ orbitals at site j defined by 

V322_^2 = \{V{rj + rdz) + V{rj - rdz)}, (1) 

and 

Vx2_y2 = V{rj + rdx), (2) 
respectively Here, V{rj) is given by 



(3) 



with a point charge ZiC at site i and the position of the 
site. rrf(=0.42A) is the radius of a Mn orbital where 
its radial charge density becomes maximum |28| and z{x) 
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FIG. 2. Tc and Tjv as functions of AV. Filled and open 
circles indicate Tc and Tjv, respectively. Al/'s are calculated 
for the same compounds in Fig. 1. Note that in the region 
with large positive (negative) AV, the 3d322_^2 (3d^2_y2) 
orbital is occupied by an electron. 

is the unit vector in the z{x) axis. The Ewald method is 
used for the lattice summation. Zj's for Mn and O ions 
and a cation at A site are chosen to be 3 -I- x, —2 and 
(8 — 2a;) /3, respectively. The difference of the potentials 



V^2_, 



(4) 



represents the relative stability of the orbitals; with in- 
creasing AV, the energy level of the Sd^^^^r'^ orbital for 
an electron relatively decreases. 

T/v and Tc are plotted as functions of Ay in Fig. |^ 
where the structural data at room temperature are used. 
Broad shades are drawn by considering experimental er- 
rors. It is clearly shown that both Tc and T/v are 
correlated with AV; T^ increases with decreasing AV 
and there is an optimal value of AV{^ 0.08 eV) for 
Tc- We estimate the strength of the correlation be- 
tween T/v and AV by using the correlation coefficient: 
r = T^Hi {Tm ~T^){AVi -~KV) / {aT^a^v) where I 
indicates a sample and N is the number of samples. 
AV (Tat) and uav (o'Tjv) are the mean value and the 
standard deviation of AV [Tm], respectively. We ob- 
tain r = —0.89 ± 0.11 for single crystal samples and 
r = —0.15 ± 0.04 for all samples including polycrys- 
tals. One might think that the Tc[n) v.s. AV curve 
in Fig. just reflects the relation between Term and x 



in La2_2a;Sri_|_22; Mn207 (LSMO) (l2[§. However, when 
we pay attention to Tc(jv)'s for samples with the same x 
(e.g., the samples C and K-R), we notice that the corre- 
lation remains. The correlation between T/v and Ay ex- 
plains that between T/v and dMn-Mn shown in Fig. 1(b), 
since AV is a decreasing function of in the re- 



gion, 3.84A < <,_M„ <3. 

In Fig. 3, we present the magnetic phase diagram at 
T = as a function of AV and x. The structural 
data at room temperature are used. Symbols connected 
by dotted lines correspond to a series of LSMO with 
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FIG. 3. The magnetic phase diagram at T = at 
function of Al/ and x. Filled, open and crossed squa 
indicate the FM, A-type AFM and PM samples, resp 
tively. AV^'s for filled and open circles are obtained fri 
the data below Tc and Tjv , respectively. Symbols conneci 
by dotted lines indicate a series of LSMO. Bold arrows sh 
5{AV) = Ay(r < Tc(iv)) - Ay(r > Tc(iv)). AF's 
calculated for the same compounds in Fig. 1. Note that 
the region with large positive (negative) AV, the 8^322. 
(3^3.2 _y2) orbital is occupied by an electron. 

X — 0.3 — 0.5 [|l2|. In addition, Al/'s calculated by usi 
the data below Tcj-tv) in LSMO are also plotted. Bold 
rows indicate the change of the Madelung potential w 
changing temperature from T > T(j(n) to T < Tc(n)'- 

5{AV) = AViT < rc(jv)) - AF(T > Tc^N^). 

We find that the magnetic structures are governed by Z 
and x; the FM (A-type AFM) phase is located in the 
gion with smaller (larger) x and moderate (smaller) A 
Let us focus on d{AV) in LSMO. S{AVys are negat 
at X = 0.3 and 0.35. The absolute value of 6{AV) grad- 
ually decreases with increasing x and S{AV) becomes a 
small positive value at a; = 0.4. Below Tc, AV seems 
to approach to the optimal value of AV^ ^ 0.08 where 
Tc becomes maximum as seen in Fig. ^ On the other 
hand, 6{AV)'s are negative at a;=0.45 and 0.48 where 
the A-type AFM structure appears. The orbital struc- 
ture and its stability in the FM phase have been studied 
by measuring the striction in Ref. The difference 

of the Mn-0 bond lengths between PM and FM states 
was reported in Ref. [Q. These experimental results are 
consistent with the present results of S{AV) in Fig. 3. 

Now we theoretically investigate a mechanism of the 
correlation between magnetic ordering and orbital struc- 
ture. We start with the following Hamiltonian |27| , ^ , p0{ : 

n = Ht + Hj + Hh + nAF + n,. (6) 

Instead of the bilayered structure, the simple tetragonal 
lattice consisting of Mn ions is considered. In this model, 
the magnetic structure with FM and AFM alignments 
perpendicular and parallel to the c axis, respectively, cor- 
responds to the A-type AFM structure . In each Mn 
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FIG. 4. The calculated magnetic phase diagram a,t T — 
as a function of A and x. Note that in the region with large 
positive (negative) A, the 3d^^2_^2 {3d^2_y2) orbital is occu- 
pied by an electron. 

ion, the two eg orbitals are introduced and the t2g elec- 
trons are treated as a localized spin with S = 3/2. The 
first two terms in Eq. (^) correspond to the so-called t- 
and J-terms in the tJ model, respectively, with the two 
Cg orbitals under strong electron correlation. The third 
and fourth terms describe the Hund coupling between 
Cg and t2g spins and the AFM SE interaction between 
t2g spins, respectively. The splitting of the energy levels 
between 3^322 _^2 and 3^3.2 _j,2 orbitals is represented by 
the last term: Ti^ = — A Tiz where the pseudospin op- 
erator is given by = ^J2-yY<y ^hcr^^'y''^'^^''^ '^^^^ 
being the annihilation operator of an eg electron at site 
i with spin a and orbital 7. The +(— ) eigenstate of Tiz 
corresponds to the state where the 3^322 _^2 (3dx2^y2) 
orbital is occupied by an electron. The anisotropies of 
the hopping integral and the SE interactions due to the 
layered structure are considered. The explicit expres- 
sion and derivation of the Hamiltonian are presented in 
Refs. 13 and 

The calculated magnetic phase diagram at T = is 
presented in Fig. 4 ]30| , |3l[| where the mean field approx- 
imation is adopted. We note that the phase diagram 
derived in this approximation explains that in the cu- 
bic manganites js^. The characteristic features shown 
in Fig. 3 are well reproduced by the present theory; the 
A-type AFM phase appears in the region with higher x 
and smaller A than that of the FM one. The range of the 
horizontal axis in Fig. 4 is larger than that in Fig. 3 by 
about 0.25. This discrepancy may be attributed to the 
neglect of the orbital fluctuation However, the char- 
acteristics of the phase diagram are insensitive to the pa- 
rameters in the model. In the FM (A-type AFM) phase, 
the orbitals are uniformly aligned with < 9 < 0.727r 
(0.727r < < tt) where 9 describes the orbital state as 
\9) —cos{^)\3d3z^-r^)— sin{^)\3d,j;2_y2) . The present re- 
sults suggest that a dimensionality of the FM interaction 
is controlled by the orbital structure; in the A-type AFM 
phase, the FM ordering in the ab plane is caused by the 
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DE interaction, while the AFM in the c direction is by the 
AFM SE. When the 3d^2_y2 orbital is stabihzed, the DE 
interaction in the ab plane (c direction) becomes strong 
(weak) and the A-type AFM phase appears |Q . A mix- 
ing of the orbitals is essential in the FM phase where 
the FM interaction overcomes the AFM SE one in the 
three directions. We note that, in Fig. 4, the FM phase 
appears not around A = but in a region of A > 0, 
since the anisotropy in the hopping integral due to the 
layered structure stabilizes the 3d^2_y2 orbital more than 
2id^^2_j.2. It is worth to mention the change of the orbital 
structure associated with the magnetic ordering: By uti- 
lizing the mean filed approximation at finite temperature, 
we compare the orbital structures above and below the 
magnetic transition temperatures. It is found that (1) 
there is an optimal mixing of the orbitals for the FM 
state and the orbital structure tends to approach to this 
structure below Tc and (2) the ?>d^2_y2 orbital structure 
is stabilized below Tjy. The theoretical results are con- 
sistent with (5(AV^)'s shown in Fig. 3 by considering that 
the change of AV^ associated with the magnetic ordering 
is caused by that of the orbital structure. 

In summary, we examine correlation between magnetic 
ordering and orbital structure in bilayered manganites. 
A relative stability of the Cg orbitals is investigated by 
calculating the Madelung potentials in a large number of 
the bilayered compounds. We find that the A-type AFM 
structure and the 3dj.2_j,2 orbital one are stabilized co- 
operatively and there is an optimal mixing between the 
3d3z2_r2 and 3d.j;2_y2 orbitals for the FM ordering. A 
theory with the two eg orbitals under strong electron cor- 
relation explains a mechanism of the universal correlation 
between orbital and magnetism. 
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